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One  weapon  system  concept  proposed  as  a  Wide-Area  Anti  armor  Munition 
(WAAM)  is  the  Antiarmor  Cluster  Munition  (ACM).  The  ACM  consists  of  a 
submunition  dropped  to  the  battlefield  by  parachute;  during  its  descent  the 
system  behaves  like  a  damped  spherical  pendulum.  At  a  certain  preprogrammed 
height  (nominally  lm)  the  submunition  detonates,  sending  forth  discrete 
warheads  at  hyper velocities.  This  method  of  delivery  gives  rise  to  certain 
ramifications  regarding  the  probability  that  a  warhead  will  hit  a  target. 
Specifically,  the  pendular  motion  acts  to  degrade  the  probability  of  hit 
PH.  This  paper  presents  a  formal  derivation  of  P^  with  regard  to  hitting 
a  solitary  target;  the  treatment  entails  coupling  the  physics  of  the  pendulum 
to  geometrical  considerations.  It  also  presents  numerical  examples  which  show 
P^  is  quite  sensitive  to  the  (ensemble-averaged)  maximum  amplitude  of  the 
pendular  motion. 
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I.  INTRODUCTION 


One  of  the  areas  receiving  widespread  attention  in  the  Air  Force  today  is 
the  development  of  Wide-Area  Anti  armor  Munitions  (WAAMs).  These  are  typically 
clustered  or  dispenser  carried  munitions  with  large  ground-plane  kill  areas 
and  are  to  be  used  in  conventional  warfare.  Several  weapon  system  concepts 
have  been  proposed  for  WAAM  applications,  one  of  which  is  the  so-called  Anti¬ 
armor  Cluster  Munition  (ACM).  This  paper  presents  a  formal  calculation  of  the 
probability  of  hit  of  the  ACM. 

The  ACM  concept  involves  the  delivery  of  clustered  bomb  units  (CBUs)  to 
the  battlefield  via  parachute.  Each  CBU  is  to  consist  of  a  discrete  number  of 
antiarmor  warheads,  or  "slugs",  to  be  explosively  dispensed  at  hypervelocity 
once  the  CBU  reaches  a  preprogrammed  height  (nominally  lm).  The  exact  con¬ 
figuration  of  the  CBU  has  not  yet  been  firmly  established;  for  concreteness, 
we  assume  it  is  shaped  like  an  equilateral  triangle  and  carries  three  slugs. 
Figure  1  portrays  the  conceptual  ACM  configuration  under  study. 

Once  the  parachute  is  deployed,  the  system  behaves  like  an  oscillator  --  a 
pendulum.  Thus,  at  the  point  of  firing,  the  system  generally  will  not  be 
oriented  vertically.  For  this  reason  the  slugs  may  fly  over  or  impact  the 
ground  short  of  targets,  and  therefore  the  pendular  motion  acts  to  degrade  the 
probability  of  hit  PH.  As  we  will  point  out,  the  behavior  of  PH  is  quite 
sensitive  to  the  orientation  of  the  CBU. 


—X  - 


Figure  1.  Conceptual  Make-up  of  the  ACM 

Our  goal  is  to  formally  derive  P^.  To  make  the  calculations  tractable, 
we  set  forth  the  following  assumptions: 

1.  The  system  always  damps  to  simple  harmonic  motion  (SHM)  before  firinq 
takes  place;  however,  the  damping  constant  itself  is  small.  This  means 
damping  over  times  as  small  as  one  period  of  the  oscillation  is  negligi¬ 
ble,  but  cumulative  damping  over  times  as  large  as  the  time  of  descent  is 
substantial. 
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2.  The  (three-faced)  submunition  always  fires  from  height  h  attempting  to 
hit  a  single  target  randomly  located  on  a  circle  of  radius  R  centered  on 
the  submunition. 

3.  The  slugs  travel  in  approximately  a  straight  line  over  their  effective 
range. 

4.  The  target  is  situated  on  flat,  open  terrain. 

5.  The  target  is  a  sphere  of  radius  h,  and  therefore  its  projection  with 
respect  to  the  submunition  is  always  the  same. 

Assumptions  2,  4,  and  5  describe  a  somewhat  contrived  situation;  however,  they 
do  not  impact  on  the  general  properties  of  P^. 

As  presented,  the  problem  is  one  of  applied  classical  mechanics,  and  its 
solution  evolves  from  the  physics  of  the  pendulum  coupled  to  geometric  consi¬ 
derations.  The  physics  of  the  pendulum  governs  the  equations  of  motion  of  the 
weapon  orientation  0,  and  these,  in  turn,  tell  how  the  orientations  are  dis¬ 
tributed  over  a  complete  oscillation;  that  is,  they  determine  the  distribution 
function  f(0)da  The  geometry  determines  the  probability  of  hit 
pH  (R,  h,  e,  $)  given  the  weapon  orientation  0  and  the  target  location  <p 
which  obeys  a  random  distribution  function  g(*)d$.  Once  pH  (R,  h,  e,  $)  is 
known,  the  calculation  becomes  purely  mathematical;  the  solution  arises  via  an 
integration  over  the  two  distribution  functions  and  becomes  a  function  of  the 
maximum  amplitude  e  .*  Finally,  an  Integration  over  a  cluster  of  ACMs 
(which  serves  as  an  ensemble),  each  member  of  which  is  deployed  using  only  one 
technique  (for  example,  deployed  from  the  same  altitude  by  the  same  type  of 
munitions  dispenser),  completes  the  calculation.  In  sections  II  and  III  these 
steps  are  developed  formally.  Section  II  treats  the  physics  of  the  pendulum, 
and  section  III  presents  the  major  aspects  of  the  calculation.  Examples  of 
the  numerical  results  are  given  in  section  IV,  as  are  the  general  conclu¬ 
sions.  Detailed  derivations  are  reserved  for  the  appendices. 

1.  pH  also  Is  generally  a  function  of  the  minimum  amplitude  s0,  but  we 
will  explicitly  treat  only  cases  for  which  0Q»o  or  0o*0m. 
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II.  PHYSICS  OF  THE  PENDULUM 


To  calculate  PH  (R ,  h,  0,  $)  we  need  only  be  concerned  with  the  state 
(orientation)  of  the  ACM  at  the  time  of  firing.  In  averaging  the  orienta¬ 
tion  9  over  its  distribution  function  f(e)d0,  an  integration  over  a  single 
period  is  required.  At  the  time  of  firing,  the  ACM  is  assumed  to  be  executing 
SHM  with  negligible  damping  over  a  single  period.  Hence,  consistent  with  our 
assumptions,  we  need  only  to  treat  the  ACM  formally  as  an  undamped  pendulum. 
However,  the  oscillations  will  generally  not  be  confined  to  a  plane;  rather, 
they  will  generally  be  "spherical".  Thus,  the  subject  of  this  section  is  the 
undamped  spherical  pendulum  executing  SHM,  and  the  goal  is  to  derive  the  cor¬ 
responding  distribution  function  f(0)d0. 

We  begin  by  deriving  the  equations  of  motion  using  the  'well-known  Lagran- 
gian  formalism  of  classical  mechanics.*  The  spherical  pendulum  is  shown  in 
figure  2;  the  submunition  of  mass  m  is  always  confined  to  a  sphere  of  radius  z 
(the  length  of  the  parachute)  and  is  subject  to  the  gravitational  accelera¬ 
tion  g.  At  time  t,  the  system  is  in  the  state  defined  by  the  amplitude  0(t) 
and  the  azimuth  y(t).  The  kinetic  energy  T  and  potential  energy  V  then  follow 
immediately  by  inspection  of  the  figure: 


(la) 


V  3  -  mgi.  cosQ 


(lb) 


The  Lagrangian  L  is  just  taT-V,  or 


L 


+  mgz  cose 


(2) 


I.  The  equations  of  motion  are  derived  in  several  textbooks,  such  as  Slater, 

J. C.,  and  Frank,  N.H.,  Mechanics,  pp.  79-86,  McGraw-Hill  Book  Company,  Inc., 
New  York,  1947. 
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Figure  2.  The  Spherical  Pendulum 


From  the  Lagrangian  follow  the  equations  of  motion;  they  correspond  to 


d  r 

3L 

3L 

-  s  n 

dt  [_3(de/dt) 

30 

-  -  u 

d  f  3L 

3  L 

-  s  0 

dt  [  3(dy/dt)  _ 

3y 

-  —  u 

(  jO  / 

By  substitution,  we  find 

- sin0  cos0  (it ■ 

dt 

x2 

■)  +  2.  si  n  9  =  o 

/  l 

(4a) 

£-  -  , 

(4b) 

where  C  is  an  integral  of  the  motion.  The  first  equation  (4)  simply  says 
energy  is  conserved,  and  the  second  says  angular  momentum  is  conserved. 

Equations  (4)  are  "exact",  but  now  we  incorporate  the  SHM  assumption  which 
says  9«1  so  that  sine  a  e  and  cose  *  1.  Then  the  equations  reduce  to 

d^Q 
dt2  ‘ 

? 

*  .  a_ 

i 

9  =  0 

(5a) 

g2  dy 

3  IT 

*  C 

(5b) 

and  upon  inserting  equation  (5b)  into  equation  (5a),  a  single  equation 
motion  for  e  follows: 

of 

d2e 

+ 

dt^ 

*-• 

c2 

'  7" 

*  0 

(6) 

The  solution  of  this  nonlinear  differential  equation  cannot  be  written  in 
terms  of  elementary  functions;  however,  it  becomes  greatly  simplified  for  two 
special  cases  -  (1)  y  *  constant,  and  (2)  e  *  constant.  For  case  (11, 
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equations  (5)  reduce  to 


■ — -f  +  -3—  9*0  0(t)  *  e  sin  uit  (7) 

dt2  4  m 

where  0m  is  the  maximum  amplitude  and  o>  -■yJg/T  is  the  circular  frequency. 
This  is  just  the  equation  of  motion  of  a  simple,  i.e.,  planar,  pendulum;  the 
orbit  in  figure  2  is  then  an  arc.  For  case  (2),  equations  (5)  reduce  to 

3  0  y(t)  3  ait  (3) 

This  is  the  equation  of  motion  of  a  "conical"  pendulum  which  executes  a  cir¬ 
cular  orbit  on  the  sphere  in  figure  2. 

Because  the  probability  of  hit  PH  must  increase  as  0  decreases,  cases 
(1)  and  (2)  in  fact  bound  the  calculations  in  that  is  maximized  for  case 
(1)  and  minimized  for  case  (2).  Thus,  we  will  focus  our  attention  on  these 
two  cases  alone.  The  derivation  of  the  distribution  function  fc(0)d8  for 
the  perturbed  conical  pendulum  is  offered  in  appendix  1,  and  it  can  be  folded 
into  the  discussion  using  the  same  techniques  presented  in  section  III. 

We  now  turn  to  developing  the  distribution  functions  for  the  simple  and 
conical  pendulums.  Equation  (7)  tells  us  how  to  form  the  distribution  func¬ 
tion  f^(e)de  for  the  simple  pendulum.  The  variable  which  determines  the 
state  of  the  system  is  the  time  t,  and  the  distribution  function  for  t  nor¬ 
malized  over  a  full  period  x  3  2ir/u  is  simply 

f (t)dt  3  (9) 

But  according  to  equation  (7), 

t  3  -  arcsin  /  —  'j  (10) 

u  \®m  / 
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so  that 


«■*  — 


vl-9 

Hence,  equation  (9)  becomes 


(ID 


fs(e)de 


(12) 


For  the  conical  pendulum,  the  distribution  function  f^ejde  is  simply  a  delta 

function  because  9 is  a  fixed  constant: 

m 


fc(e)de 


dr 6(Q  *9m  )d0 


(13) 


Both  f^(9)  and  ■f(-( Q)  are  normalized  over  a  full  period: 


fs(0)de 


3 m 

4/fc 


.(©)d@ 


1 


(14) 
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III.  THE  PROBABILITY  OF  HIT 


The  probability  of  hit  PH  (R,  h,  a,  <t>)  given  the  orientation  9  of  the  ACM 
and  the  azimuth  $  of  the  spherical  target  arises  solely  from  geometrical  con¬ 
siderations.  The  calculation  of  the  expected  probability  of  hit 
PH  (R»  h,  <0m>)  involves  integrations  over  the  random  distribution 
function  g($)di>  and  the  (not  random!)  distribution  function  f(9)de  derived  in 
the  last  section.  Thus,  we  begin  this  section  with  a  discussion  of  the 
geometry  and  then  proceed  with  the  calculations. 

A.  GEOMETRY  OF  THE  ACM  ATTACK 

When  the  submunition  fires,  its  plane  of  fire  generally  is  not  parallel  to 
the  plane  of  the  ground.  This  means  that  in  certain  cases,  i.e.,  for  certain 
R  and  9,  the  slug  may  fly  over  the  target  or  may  impact  the  ground  short  of 
the  target.  This  situation  Is  portrayed  in  figure  3  which  shows  the  weapon's 
plane  of  fire  is  inclined  by  the  angle  0  to  the  "equator"  which  runs  parallel 
to  the  ground  and  through  the  submunition  firing  at  height  h.  The  figure 
makes  clear  that  the  target  dimensions  give  rise  to  a  "lethal"  region  bounded 
by  the  "critical  azimuths"  ±<f>CR.  (There  are  actually  two  lethal  regions 
which  are  diametrically  opposed.)  If  0< a/2,  the  angular  extent  of  the  target, 
then  the  target  certainly  lies  within  the  plane  of  fire,  i.e.,  $CR  *  *72.  The 
probability  that  a  slug  hits  a  target  outside  the  lethal  region  is  zero. 

As  figure  3  implies,  the  calculations  involve  spherical  geometry  because, 
from  the  submunition's  viewpoint,  it  is  the  target's  projection  on  the  sphere 
of  radius  R  (corresponding  to  the  range)  that  is  being  shot  at.  The  angle 
subtended  by  a  spherical  target  of  radius  h  at  range  R  is  illustrated  in 
figure  4  from  which  we  infer 

a  ■  2  arcsin  ( jp)  (15) 
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ZENITH 


WEAPON'S  PLANE 
OF  FIRE 


GROUND  LEVa 


"LETHAL"  4ro  "  CRITICAL 
REGION  U  ANGLE 


Figure  3.  Geometry  of  the  ACM  Attack 


Figure  4.  Target  Geometry 

Obviously,  for  3»0,  the  probability  of  hitting  a  single  target  at  range  R, 
assuming  no  other  targets  or  obstacles  lie  inside  the  circle  of  radius  R,  is 
simply 

ph  (r,  h,  3*o)  » |f  *4"  arcsin(ir)  (16) 

Here,  the  factor  3  arises  because  we  have  assumed  the  submunition  simul¬ 
taneously  fires  three  slugs  at  azimuths  separated  by  120°.  Clearly,  Ph 
drops  from  its  maximum  value,  unity,  at  R*h/sin(ir/3). 

The  geometry  of  the  lethal  region  as  viewed  from  the  submunition  is  lifted 
out  of  figure  3  and  shown  In  figure  5.  If  the  target  lies  inside  ±^,  then 
there  Is  a  non-zero  chance  for  a  hit,  but  if  the  target  lies  outside  this 
boundary,  there  is  no  chance  for  a  hit.  This  is  precisely  what  figure  5  shows. 
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F 


Figure  5.  Geometry  of  the  Lethal  Region 

Vie  now  derive  from  figure  5  using  the  spherical  triangle  ABC  which 
arises  from  the  target  dimensions  and  the  orientation  of  the  ACM.  The  angle  v 
is  denoted  merely  as  an  auxiliary  angle  to  assist  in  the  derivation.  The 
spherical  triangle  relations  which  follow  from  the  figure  are: 


sini|>  =  cos  esec(a/2) 

(17a) 

cos  $CR  *  cot  Scot  <l> 

(17b) 

Combining  these  gives 

*CR  a  arcsin  [  s'in^1  ] 

(18) 

Equation  (18)  holds  for  0  >  a/2;  If  0  <  a/2,  then  *CR  *  ir/2  always.  Under 
the  SHM  assumption,  the  angle  a/2  is  always  small  when  0  >  a/2,  but  it  is  not 
necessarily  small  when  S  <  a/2.  Thus,  under  the  SHM  assumption  and  according 
to  equations  (15)  and  (18),  the  critical  angle  is  given  by 

*CR  *  ! 

ir/2  for  0  <  a/2 

arcsin  (h/R©)  for  0  > 

*  h/R 

(19) 

a/2  ■  arcsin  (h/R) 

v 
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8.  DERIVATION  OF  THE  PROBABILITY  OF  HIT 


We  know  that  if  the  target  lies  outside  the  lethal  region  of  figure  3, 
then  the  corresponding  probability  of  hit  is  zero.  Otherwise,  however,  PH 
is  non-zero,  and  we  now  seek  to  quantify  this. 


We  begin  by  examining  figure  6  which  shows  a  projection  of  the  target 
lying  inside  the  lethal  region.  The  angle  exposed  to  the  ACM’s  plane  of  fire 
is  S,  and  the  angles  X  and  are  auxiliary  angles  to  assist  in  deriving  6. 
Once  8  is  known,  then  the  probability  of  hit  is  known;  it  is  just  3e/2ir. 

We  now  derive  0  using  the  spherical  triangles  ACD  and  BCD  in  figure  6. 

The  relations  corresponding  to  triangle  ACD  are 

cos  X  3  cos  9  CSC  rfr  (20a) 

tan  1/  *  cot9  sec  <P  (20b) 

and  the  relation  corresponding  to  triangle  BCD  is 

cos  (s/2)  *  cos  (a/2)  sec  x  (20c) 


Combining  these  gives 


8 


2  arccos 


LV 


COS  (a/2) 


2  2 
1  -  sin  0  sin  4 


(21) 


Thus,  the  probability  of  hitting  a  single  target  at  range  R  given  that  the 
target  lies  at  azimuth  <p  inside  the  lethal  region,  and  given  that  the  orienta¬ 
tion  of  the  ACM  is  9,  is 


PH  (R,  h,  0,  *) 


36  3 

*  —  arccos 

2*  IT 


l-(h/RP 

2  2 
1-sin  0  sin  $ 


(22) 


The  comments  made  in  conjunction  with  equation  (16)  also  apply  here.  Note 
that  for  9*0,  equation  (22)  reduces  to  equation  (16). 
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ALL  ARE  EXPOSED  ANGLE  8 


Figure  6.  The  Exposed  Target 

Next  we  need  to  average  P^  (R,  b,  9,  $)  over  the  target  azimuth  <j>  and  the 
ACM  orientation  e;  we  begin  by  averaging  over  Because  the  orientation  of 
the  submunition  is  random,  $  is  a  random  variable,  and  its  distribution 
function  g(<fr)d$  is  therefore  random: 

g(*)d$  3  (23) 


To  compute  Ph  (R,  h,9)  requires  Integrating  equation  (22)  over  g(<t>)d$: 


Ph  (r»  h>  ®) 


[PH  (R.  h,  0, 

■</ 

SO 


$  g( <P)dg> 


<tR 

PH  (R,  h,e, 4>)  g(<o)d<(» 


(24) 
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By  substitution,  we  have 


arcsin  (h/Re) 
PH  (R,  h,  0)  *  -4y-  J  arccos 


l-(h/R) 


r 


5  .2 

1  -  sin  0  sin  4,' 


dtp  (25) 


Using  the  SHM  assumption,  we  can  solve  this  integral  in  closed  form,  and  we  do 
so  fn  appendix  8.  The  solution  corresponding  to  ^  *  ff/2  is 


PH  (R.  h,  e<  h/R)  *  i- 


arccos 


L  cose 


lrV1'®?[E(Re/hl '  fV'W 


and  the  solution  corresponding  to  <t^<ir/2  is 


PH  (R,  h,  9>  h/R)  *  nr  [  “T  E(h/Re) 

it  R  L. 

-(t*  -  fe)  K<h^>] 


(26a) 


(26b) 


In  these  equations  K  and  E  are  the  complete  elliptic  integrals  of  the  first 
kind  and  second  kind,  respectively,  and  are  elementary  functions.^ 


The  goal  now  is  to  average  PH  (R,  h,  e)  over  the  distribution  functions 
f(.(e)de  and  f^(9)de  for  the  conical  and  simple  pendulums.  Because  ) 
is  a  delta  function  (see  equation  (13)),  the  integration  is  Immediate,  and  for 
the  conical  pendulum  9  is  replaced  by  sm  in  equations  (26).  What  is  left, 
then,  is  to  evaluate  the  following  integral  with  regard  to  the  simple  pendulum: 

(o)ds  (27) 


P„  (R,  h. 


9m> 


L 


m 


P^  (R,  h,  0)  f  ^ 


1.  Abramowltz,  M.,  and  Stequn,  I. A.,  Handbook  of  Mathematical  Functions,  pp. 
589-626,  National  8ureau  of  Standards,  1964.  (Short  Title  -  Abramowitz  and 
Stegun) 
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The  integral  cannot  be  solved  in  closed  form;  as  is  shown  in  appendix  C,  the 
solution  leaves  three  integrals  to  be  evaluated  numerically.  For  *  ir/2, 
the  solution  is 


and  for  $CR<ir/2,  the  solution  is 


Pjj  (R,  h,  9  ^  >  h/R) 


<V*,> 


where 


(28a) 


(28b) 


(28c) 


(28d) 


(28e) 
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Finally,  for  completeness,  we  must  average  PH  (R,  h,  0m)  over  an 

ensemble  of  ACMs,  each  member  of  which  is  deployed  using  the  same  technique, 

to  calculate  P„  (R,  h,  <©>),  where  <0  >  is  the  average  maximum  amplitude 

to  be  determined  observationally.  For  a  less  restricted  treatment  in  which 

the  SHM  assumption  is  not  made,  an  integration  over  a  Gaussian  distribution  of 

tan^  would  now  be  appropriate,  and  the  result  would  be  PH  (R,  h,  <9^  ,  a) 

where  a  is  the  standard  deviation  of  the  <0m>  observations.  However,  we  have 

assumed  the  ACM  always  executes  SHM  when  it  fires,  and  in  doing  so  we  have 

essentially  assumed  3  <<  <  ®m  ><  Thus,  we  have  a  leptokurtic  (strongly  peaked) 

distribution  function  for  centered  on  <0  >,  and  consistent  with  the  SHM 

m  m 

assumption  we  can  make  an  approximation  and  simply  replace  by  <©m>  in  the 
previous  equations. 


IV.  RESULTS  AND  CONCLUSIONS 


We  now  discuss  the  sensitivity  of  P^  to  the  average  maximum  amplitude 
<em>,  the  target  dimension  h,  and  the  mode  of  oscillation  (simple  versus 
conical).  Figures  7,  8,  and  9  were  developed  using  available  polynomial 
expansions  for  the  complete  elliptic  integrals*  and  straightforward  Riemann 
summation  to  perform  the  integrals  in  equations  (28).  As  examples  we  have 
selected  <0  >  *  1°,  4°,  7°,  and  10°  and  have  arbitrarily  restricted 
our  curves  to  ranges  R  £  30m. 


Figure  7  shows  four  P^(R )  versus  R  curves,  one  for  each<em>, 

plotted  on  a  log-linear  scale.  The  mode  of  oscillation  is  that  of  a  simple 

pendulum  and  h*lm.  The  curves  begin  to  drop  from  unity  at  R*h/sin(ir/3)  and  as 

R  increases,  the  spread  between  the  four  curves  also  increases  and  quickly 

becomes  substantial.  For  example,  at  R*30m,  the<^T)>*  1°  and<^>  ■  4° 

curves  differ  by  a  factor  of  1.8,  and  the  1°  and  10°  curves  differ  by  a 

factor  of  3.6.  The  slight  "knees"  circled  in  the  4°,  7°,  and  10°  curves 

mark  the  points  at  which  <0fJ1>  *  h/R;  P^( R)  falls  off  slightly  more  rapidly 

once  <9>  exceeds  h/R. 
m 


In  figure  8,  the  sensitivity  of  PH  to  the  mode  of  oscillation  is 
illustrated.  Both  P^(R)  and  pj^(R)  versus  R  are  plotted  for<em>*  10° 
and  h*lm.  Again  there  is  a  rapidly  growing  spread  as  R  increases.  At  R*30m 
the  curves  differ  by  a  factor  of  2.8.  Not  shown  are  the  <em>  s  1°,  4°, 
and  7°  examples;  for  these  cases  the  curves  differ  by  factors  of  1.04,  2.1, 
and  2.4,  respectively. 


1.  Abamowitz  and  Stegun,  equation  (17.3.34),  p.  591,  and  equation  (17.3.36), 
p.  592. 

2.  At  small  R  the  submunition  can  conceivably  collide  with  the  target  before 
firing;  this  fact  is  not  accounted  for  in  the  numerical  examples. 


Figure  7.  P{j(R)  Versus  R  -  The  Simple  Pendulum  Mode 
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RANGE  (m) 


Figure  8.  Simple  Pendulum  Mode  Versus  Conical  Pendulum  Mode 

s 

Finally,  figure  9  shows  the  sensitivity  of  PH  to  h.  Examples  are 
plotted  which  correspond  to  o  >*  1°  and  10°,  and  h*lm  and  1.5m.  The 
spreading  between  the  <0m>  ■  1°  and  <©m>  s  10°  curves  also  increases 
rapidly  for  h*1.5m,  but  at  R*30m,  the  curves  differ  by  a  factor  of  2.7.  Thi 
is  less  than  the  factor  of  3.6  for  h»lm  but  is  nevertheless  quite  subtantial 


ph«r) 


Figure  9.  Pj^(R)  Versus  R  -  Sensitivity  to  Target  Dimensions 

The  discussion  thus  far  leads  pointedly  to  a  single  conclusion:  Against 
targets  that  the  ACM  would  typically  encounter  on  the  battlefield,  the  proba¬ 
bility  of  hit  is  sensitive  to  the  average  maximum  amplitude.  The  sensitivity 
becomes  more  pronounced  with  larger  om>»  especially  with  regard  to  a  target 
at  a  large  range.  The  sensitivity  is  further  accentuated  if  the  mode  of 
oscillation  deviates  from  that  of  a  simple  pendulum. 
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At  this  time  we  point  out  that  we  have  treated  only  P^  and  not  the  pro¬ 
bability  of  kill  PK,  for  PK  is  a  product  of  two  factors  —  (1)  the  proba¬ 
bility  of  kill  given  a  hit  P(K|H),  which  depends  on  target  characteristics, 
and  (2)  the  probability  of  hit  P^: 

Pk»P(K|H)Ph  (29) 

An  analysis  of  ACM  performance  in  a  realistic  scenario  must  therefore  also 
include  P(K|H).  Nevertheless,  we  can  still  conclude  that  a  highly  effective 
and  reliable  stablizing  mechanism,  or  possibly  a  homing  device,  is  required  in 
the  ACM  design  to  avoid  the  risk  of  severe  degradation  of  P^  and  therefore 


APPENOIX  A 


DISTRIBUTION  FUNCTION  OF  THE  PERTURBED  CONICAL  PENDULUM 


Although  the  equation  of  motion  (6)  of  a  spherical  pendulum  executing  SHM 
cannot  be  solved  in  closed  form,  we  nevertheless  can  investigate  a  certain 
restricted  case  analytically;  namely,  the  perturbed  conical  pendulum.  For 
this  case  there  are  both  a  maximum  amplitude  0  and  a  minimum  amplitude  9 

m  0 

such  that 


0 

m 


0 


(Al) 


By  use  of  condition  (Al)  and  a  Taylor  series  expansion,  we  will  convert  the 
equation  of  motion  (6)  into  a  differential  equation  for  0  that  can  easily  be 
solved.  We  will  then  convert  the  solution  into  the  distribution  function 

'V, 

f(.(9)d9  of  the  perturbed  conical  pendulum. 


The  equation  of  motion  of  the  spherical  pendulum  executing  SHM  is 


*  JQ  .  _cL  *  0  (A2) 

dt  9J 


Consistent  with  condition  (Al),  we  transfer  the  time  dependence  of  0  over  to  a 

perturbation  variable  n«0  : 

m 

©(t)  *  em  -  n(t)  (A3) 

Then  a  Taylor  expansion  to  first  order  in  n  gives 


Substitution  of  equations  (A3)  and  (A4)  Into  equation  (A2)  gives  a 
second-order  linear  differential  equation  for  n(t): 


where 


d2n 


dt 


*“!"•  °fa«  ■  0 


(A5) 


w2  +  3 

V 


(A5a) 


=  “2 


C 


(A5b) 


To  solve  equation  (A5),  we  demand  that  the  following  boundary  conditions  hold: 


n ( t*o)  *  0 


(A6a) 


JL^SI  ,0 


(A6b) 


The  physically  meaningful  solution  satisfying  these  conditions  is 

n(t)  »  [l-  |cos  (a  jt)|]  (A7) 

Thus,  the  perturbation  varies  periodically  with  a  period  t»  2 Tr/rz ^ . 

There  Is  still  one  unknown  constant,  C,  to  evaluate,  and  therefore  we  need 
one  more  boundary  condition.  Consistent  with  the  boundary  conditions  (6),  we 
demand  that 


n(t  «  t/4)  »  0-0 

m  o 


(A8) 
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ft 


This  gives  from  equation  (A7) 


*  "  0« 
m  o 


and  by  substitution  of  equations  (A5),  we  have 


(U2  -  /-V)  0. 


<Ji2  +  3/— (0  -  0  ) 

l  02  I  m  0 

\  m  /  • 


from  which 


C  3  9m  “\f  ♦  3?0m-eo] 


(AlO) 


2  2 

The  expressions  for  Qj  and  ^  now  reduce  to 


4V 

\  *  3(V9o> 


(Alla) 


.2  _  4(V9o)m 

2  ^  ♦  3(0m-0o) 


(Allb) 


and  the  expression  for  e(t)  reduces  to 


e(t)  »  s0  +  (0m  -  e0)  I  cos^t)  | 


(A12) 


If  we  rewrite  equation  (A12)  in  the  form 


0{t)  ■  0Q  +  (0m  -  @0)  i  sinfajt  +  -J-)  | 


(A13) 
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then  it  is  straightforward  to  show 


it  ■  h  (A14) 

V  (vv2  •  (9-9°>2 

The  distribution  function  f^(0)d©  then  follows  in  the  same  way  as  equation 
(12)  of  section  II: 


?c(e)de 


2ir 


_ d0 _ 

V  (6»  -  9°>  •  ,S  -8o>2 


(A15) 


This  distribution  function  can  be  folded  into  the  calculations  in  section  III 
directly,  except  now  two  observed  quantities  are  involved  —  <0m>  and 

<0Q>. 


APPENDIX  B 

THE  CALCULATION  OF  PH  (R,  h,  0) 


This  appendix  details  the  calculation  that  carries  equation  (25) 
equations  (26).  Equation  (25)  is 


arcs  in  (h/R 0) 


PH  (R,  h,  0  )  *  — j  I  arccos 


f 

J  jl-sin^Q  sin 


d<$> 


To  simplify  the  notation,  we  rewrite  this  in  terms  of  a  and  <p. 


CR‘ 


PH  (a’  9*  *CR}  *  7  J 


'CR 

arccos 


cos  (g/2) 
1-sin^a  sin2$  J 


d<|) 


An  Integration  by  parts  gives 


f 

»  _ 

(a,  0,  $cr)  *  ~r 

arccos 

cos  (a/2) 

L  .2.2 
Wl-sm  esm  <j>CR 

♦  cos  (a/2) 


SCR 

7. 


$sin$cos$  sin  0d<l> 

- 5 - 

.‘a 


(1-Sin  Q  sin  $)  ^  sin2(a/2)  -  sin20si n% 


We  now  use  the  $HM  assumption  to  reduce  the  remaining  Integral  which 
"I": 


I  (a,  0,  $CR) 


J  J 


V  sin2(a/2)  -  02  sin24> 


over  into 


(Bl) 


(32) 


(B3) 


we  call 


(B4) 
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and  an  integration  by  parts  gives 


To  reduce  the  remaining  Integral  we  again  change  variables  and  introduce  a  new 
parameter 

It  *  Q  esc  (a/2)  ■  RS/h  (B7) 


which  results  in 
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The  remaining  integral  is  the  elliptic  integral  of  the  second  kind1,  and  so 
we  have  the  result 


PH  (a’ 0*  W  SA 


?CR 


arccos 


.7 


cos (a/2) 


2  2 
1-sin  0sin  a 


crj 


sin« 


E( *CR  ’  ‘  *CR^(*"k^  S1'n"  ^CR 


(B9) 


For  <(>cr  3  *72,  this  reduces  to 

♦CR 


PH  (a,  0,  ?CR  *  ir/2)  3  T 


arccos 


E 


CQS( g/2) 
cose 


-VTTTl 


(810) 


where  E(k)  is  the  complete  elliptic  integral  of  the  second  kind.  Equation 
(26a)  follows  immediately  upon  substituting  the  expression  for  k  (equation 
(87))  and  for  a  (equation  (15))  into  equation  (BIO).  Note  that  for  s=0, 
equation  (BIO)  reduces  to  equation  (16). 


To  recover  equation  (26b),  let  us  first  look  at  EC**^,  k): 


E(* 


CR 


,s,n*CR - . 

•k>5/ 


d  y 


(BU) 


1.  Gradshteyn,  I.S.,  and  Ryzhlk,  I.W.,  Table  of  Integrals,  Series,  and 
Products,  p.  905,  Academic  Press,  Inc.,  New  York,  196b.  (Short  Title  - 
Gradshteyn  and  Ryzhlk)  Note  that  Gradshteyn  and  Ryzhik's  k«  is  Abramowitz 
and  Stegun's  m. 

2.  Gradshteyn  and  Ryzhlk,  p.  905. 


According  to  equation  (19) 


sin$££  *  h/RQ.  *  1/k 
Thus,  equation  (811)  becomes 


(B12) 


E(l/k,  k)  * 


(B13) 


A  change  of  variables  gives 


E(l/k,  k) 


This  integral  can  be  solved  in  closed  form;  the  solution  is^; 


(B14) 


E(l/k,  k) 


k  E( 1/k )  - 


K ( 1/k) 


(B15) 


2 

where  K  is  the  complete  elliptic  integral  of  the  first  kind.  If  we  now 
substitute  equations  (18)  and  (B15)  into  equation  (B9),  we  find  (remembering 
sin0*0) 


(a »  0»#£h<it/2) 


3  sina 


■  k  E(  1/k )  -  K(l/k) 


(B16) 


Equation  (26b)  now  follows  immediately  upon  substituting  the  expressions  for  k 
and  o. 


1.  Gradshteyn  and  Ryzhlk,  equation  (3.169.9),  p.  276. 

2.  Gradshteyn  and  Ryzhik,  p.  905. 
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APPENDIX  C 


THE  CALCULATION  OF  PH  (R,  h,  e^) 

This  appendix  details  the  calculations  that  carry  equation  (27)  over  into 
equations  (28).  We  begin  by  deriving  equation  (28a).  Substituting  (12)  and 
(26a)  into  equation  (27)  gives 


An  integration  by  parts  of  the  first  integral  gives 


Using  the  SHM  approximation,  this  becomes 
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Another  Integration  by  parts  gives 


Equation  (28a)  follows  Immediately  upon  substituting  this  into  equation  (Cl). 

We  now  turn  to  equation  (28b).  For  the  case  0m  >a/2  two  terms  arise. 

The  first  corresponds  to  replacing  3m  by  a/2  in  equation  (Cl),  and  the 

second  corresponds  to  an  Integration  of  equation  (26b)  over  f^(0)dQ  from 

«/2  to  0„.  Thus,  we  have 
m 


A  change  of  variables  to  x  *  arcsin  (©/0m)  results  in  equation  (28b).  The 
second  integral  can  be  reduced  to  two  integrals  involving  only  E  and  not  K  in 
their  kernels,  but  these  still  require  numerical  integration. 


Note  that  for  the  case  9m  ■  h/R  equations  (28a)  and  (28b)  match: 
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The  integral  has  a  closed  form  solution,!  and  the  result  is 


1.  Gradshteyn  and  Ry2hik,  equation  (6.151),  p.  638. 
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LIST  OF  ACRONYMS  AND  SYMBOLS 

ACM 

Antiarmor  Cruise  Munition 

C 

Integral  of  the  Spherical  Pendulum  Motion  (see  equation  (4b)) 

E 

Complete  Elliptic  Integral  of  the  Second  Kind 

f 

Distribution  Function 

fs(0) 

Distribution  Function  of  the  Orientation  0  of  the  Simple  Pendulum 

fc(0) 

Distribution  Function  of  the  Orientation  q  of  the  Conical  Pendulum 

V®> 

Distribution  Function  of  the  Orientation  e  of  the  Perturbed  Conical 

Pendulum  • 

g 

Gravitational  Acceleration 

h 

Radius  of  Spherical  Target 

i 

Auxiliary  Integral  (see  equation  (B4)) 

!-i  •  I2*  I3 

Integrals  Comprising  P^  (see  equations  (28)) 

Auxiliary  Integral  (see  equation  (C2)) 

k 

Auxiliary  Parameter  (see  equation  (B7)) 

K 

Complete  Elliptic  Integral  of  the  First  Kind 

l 

Length  of  Parachute 

L 

Lagrangian 

m 

Mass  of  Submunition 

Probability  of  Hit 

PS 

Probability  of  Hit  of  the  Conical  Pendulum  Model 

* 

Probability  of  Hit  of  the  Simple  Pendulum  Model 

PK 

Probability  of  Kill 

P  ( K | H) 

Probability  of  Kill  Given  a  Hit 

R 

Range 

SHM 

Simple  Harmonic  Motion 

t 

Time 

T 

Kinetic  Energy 

7 

Velocity 

V 

Potential  Energy 

WAAM 

Wide-Area  Antiarmor  Munition 

a 

Maximum  Angle  Subtended  by  Target  (see  equation  (15)) 

3 

Angle  Subtended  by  Exposed  Target  (see  equation  (21)) 

41 

Azimuth  of  Spherical  Pendulum 
Perturbation  Amplitude  (see  equation  (A3)) 

ACM  Orientation  (see  figure  2) 

Maximum  Amplitude  of  Orientation 
Average  Maximum  Amplitude  of  Orientation 
Minimum  Amplitude  of  Orientation 
Standard  Deviation 
Period  of  Pendulum 
Azimuth  of  Target 

Auxiliary  Angle  (see  equations  (20)) 

Auxiliary  Angle  (see  equations  (20)) 

Circular  Frequency  of  Pendulum 

Circular  Frequency  of  Perturbed  Conical  Pendulum 

Auxiliary  Frequency  (see  equation  (A5b)) 


